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Recently, the free energy of the target space mean field (TSMF) matrix model has been 
calculated in the low temperature phase, order-by-order in a low temperature expansion. 
The TSMF model is a matrix model whose discrete target space has an infinite coordination 
number, and whose free energy assumes a universal form, corresponding to baby universes 
joined into a tree. Here the free energy is summed to all orders, and expressed through 
a transcendental algebraic equation, using which we analyze the critical phenomena that 
occur in the TSMF model. There are two critical curves, at which the matter and the 
geometry become critical, and for which the critical exponents are ct = ^ and 7str = 
respectively. There is a bicritical point where the curves meet, at which 7str = 
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Matrix models — convenient representations of two-dimensional random surfaces and 
string theory — have run into a brick wall, the central charge c = 1 barrier. The models are 
exactly solvable when the target graph is a tree (that is, it has no loops) and the matter 
has c < 1, and almost completely mysterious otherwise. 

Recently it has been noticed that matrix models on target graphs that are regular and 
whose coordination number A is large are rather simple, at least in their low temperature 
phase These models have target graphs with many loops (for example, (A + l)-state 
Potts models) , and often positive infinite central charge (log^ A uncoupled g-state Potts 
models, and A/2-dimensional lattices). In particular, the leading 1/A term in the free 
energy contains contributions from all labeled trees, where each vertex of each tree is a 
pure-gravity, open baby universe. Furthermore, this leading term is universal: it is entirely 
independent of the target graph (and therefore of the matter model or embedding space) . 
This is rather reminiscent of mean-field theory on ordinary lattices, where the system 
becomes independent of the particular matter model when the (world-sheet) coordination 
number gets large. For this reason we call the present universal matrix model "target space 
mean field theory" (TSMF). 

We begin by briefly summarizing the results of the previous work |jl|] . The free energy 
of the TSMF model is a function of the cosmological constant g and of a parameter a, which 
depends on a = ^ where (3 is the inverse matter temperature of the matrix model. For 
a single matter model, for example, a = Aa; for multiple matter models, a = e^" — 1, 
where now A is a product of the number of models and the coordination number of each 
model. The free energy was expressed as a low-temperature expansion, a power series in 
a: 



where 7^ is the set of all labeled free trees with n vertices. (For graph theoretical deflnitions 
and examples, see, for example, 0.) At this stage, the reader might worry that a is not 
small if A is large; this is compensated by the fact that for large A, the critical temperature 



each tree T, V{T) is a product of a vertex factor Pn{g), for each n-degree vertex in T; for 
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so that Qc remains flnite, and (^ is valid all the way to the singularity 0. For 




Thus, for the flrst few orders. 
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The vertex factors are just moments of a hermitian one-matrix integral: 

A=0 



Pnig) 



where = D(l)e-^'^^^\ U{<p) is a non-gaussian matrix-model potential, and (p is an 

N X N hermitian matrix. In this work we will treat the cubic potential 

U{<f>) = 1<P' + -^<P' (4) 



which generates triangulated surfaces, and at the end mention other, higher-degree, poten- 
tials. The generating function Il{g, A) can be evaluated by the methods of 0], for instance; 
we will give explicit expressions for it below. We shall use the following notation for the 
partial derivatives of 11 : 

niU,A) = — — — , U2{g,X) = — ^ — , etc. (5) 

The free energy series (|1|) is the leading, A° term in the 1/ A expansion, and is a great 
simplification over the general low-temperature series The subleading terms (see for 
examples) are neither tree-like nor universal. 

In we gave closed-form expressions for each order in the low-temperature expansion 
(p. In order to analyze the critical behavior of the free energy, what is needed is a sum of 
the entire series. To calculate it, we introduce the scalar integral 

z{pi, a, e) = ^r^^ (6) 

/ ax e ■i 

J — oo 

This is a generating function for closed graphs with labeled vertices of all degrees, a degree 
i vertex receiving a factor ^^,0 and the graph with n vertices and £ lines receiving a factor 
^nQn-£_ The highest possible power of 9 will be 9^, which will multiply trees. By taking 
the logarithm and the 9 ^ oo limit, and then identifying the coefficients pi with the vertex 
factors of (^), we obtain exactly the TSMF free energy: 

F(^,a) = - lim ^log^(a,^) (7) 



^ We need to divide pi by i\, lest each line out of each vertex be individually labeled. 
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We have thus integrated out the hermitian matrices of the multi-matrix model, leaving 
behind a single scalar. In the original matrix model, each vertex represents a point on a 
surface; in the effective theory defined by (||) and (|^), a vertex represents a baby universe. 

Two comments are in order here. First, the F in (|^) is not really a sum of the series 
(P, but a function whose low temperature expansion coincides with (|l|). The free energy 
is expected to be analytic at a = 0. Therefore we can trust up until a singularity. 
Thereafter, as we shall see below, the 9 —>■ oo limit is poorly defined, as the action no 
longer has a saddle point. 

Second, it would be exceedingly pleasant if the subleading terms in 1/6' in (|^) cor- 
responded to the subleading terms in 1/A in the free energy of the matrix model. As 
we pointed out in [Q, the latter are no longer universal; but perhaps some matrix model 
could be found for which the subleading — and not just the leading — terms in the 9 and A 
expansions coincide, and for which, therefore, the free energy in (0) but without the limit 
would be the exact free energy, to all orders in 1/A. Unfortunately, such a matrix model 
is not likely to exist. This is because in the limit A ^ oo, the baby-universe vertex factors 
Pn (0) include contributions from twisted surfaces (see for an example). As soon as we 
add loops to the tree, some of the configurations with twisted vertices are suppressed, as 
they result in non-spherical surfaces; therefore the vertex factors become dependent on the 
global properties of the tree, and cannot be expressed by means of a local theory such as 
(^, (0). Perhaps this difficulty could be circumvented by taking the double-scaling limit, 
in which all genera are present. 

Now, back to the calculation. We evaluate (^) by steepest descents, of course. Our 
final expression for the free energy is 

a F{g, a) = a n(^, y{g, a)) - ^y{g, af (8) 

where the function y{g, a) must be a local extremum which maximizes the free energy: 

ani(5f, y(5f,a)) - y(5f,a) = (9) 

In the same way that F{g, a), when expanded in powers of a, generates free trees, y{g, a) 
generates planted trees, that is trees with one chosen, degree-one vertex. We will call (^), 
our "equation of motion," the tree equation. The low temperature expansion of y begins: 

y{9, a) = Pia + PiP2a^ + {pivl + «^ ^ (10) 
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Equipped with the convenient form (^) of the free energy, let us examine its "physically 
relevant" critical regime(s). Let us first recall the general features of the critical behavior 
of the c < 1 models. There, the surface does not yet break up into baby universes, and 
the temperature is coupled to the order parameter of the matter. Given a generic value 
of the temperature, there is a critical value of the cosmological constant which makes the 
area of the surface diverge; generically, 7str = There is, in addition, one bicritical 

point, at which the matter also undergoes an ordering phase transition, which modifies the 
geometric transition, so that 7str is between — i and 0. 

In the case of the TSMF model, we expect somewhat different behavior [Q. The cos- 
mological constant g is still coupled to the area of the surface, but as it only appears in the 
vertex factors, it is coupled to the area of the baby universes, rather than the entire surface. 
The (transformed) matter temperature a is no longer coupled to the "magnetization" of 
the matter; instead, as can be seen from (|I|), it is coupled to the size (number of nodes) of 
the tree. In this sense it becomes a second-quantized cosmological constant, coupled to the 
number of baby universes. We still expect a phase transition, no longer a matter-ordering, 
but a tree-growing phase transition, the point where the number of baby universes di- 
verges. Thus, there will be two critical curves in the {g, a) plane: the "T-curve," where the 
number of nodes (baby universes) in the tree diverges; and the "G-curve," the geometrical 
phase transitions, where the area of each baby universe diverges. If the curves cross there 
will be a bicritical point, at which both the number and the area of the baby universes 
diverge. 

We define the exponents as closely as possible to the usual conventions. Let d be the 
perpendicular distance in the {g, a) plane between some point and the T-curve; then as 

F(£?, a) ^ analytic + 5^"" + ■■ ■ (11) 

We have chosen the symbol a because the T-curve is a remnant of the matter-ordering 
transition. Similarly for the G-curve: let 5 be the deviation from it, and define 

F{g, a) ^ analytic + J^-T^tr ^ . . . (22) 
Of course, at the bicritical point the two exponents coincide. 



4 



The first thing to do is to calculate the generating function n((7, A), which is easily 
done by shifting the matrix (p ^ (p + xl, in order to eliminate the linear term; this gives 

U(g, A) = -A(x + A) - - A^) -gx^-^ log(l - 12gX) 

-n(r(<7,A)) 

-1 + VI - 12gX 

^ = a 

6g 

T{g,X) = {l-12gX)-^/^g 

and n(^) = -n(^,0).i An explicit expression for Il{g) can be found in [||, and is repro- 
duced in the Appendix. The geometric singularity in Il{g) occurs at g = go = 1 / a/ 108a/3 ~ 
0.0732. 

The G-curve is the easiest to calculate. The transition occurs when r{g, z/q) = go; in 
other words, 

We then solve the tree equation to obtain ao{g). Let us analyze the critical region by 
expanding the tree equation; putting Aa = ao — a, Ay = yo — y, we get 

uig) + v{g)Ay + w{g)Ay^'^ + . . . = ^ _ ^ + ^Aa + . . • (15) 

ao ao ap 



where the functions u, v, and w can be determined from (|T3D after a little algebra: see 
Appendix. We solve for the critical points by canceling the analytic parts of ([T5|) . Canceling 
the constant term 

Go = yo/u (16) 

gives an explicit form for ao = 9 / 9o)- 

1 - 7^/3 

" 1 - 272/3 + 74/3 + 36(^2^2/3 + 15552 . 31/4 (3 + y3)-3^3^2/3 '^^^^ 
If we also cancel the linear term, that is if 

ao{9)v{g) = -1 (18) 



2 Note the sign convention: all coefficients in the expansion of Il{g, A) about (0, 0) are positive. 
Therefore our free energy at a = will have a negative sign with respect to the function E^^\g) 



m 



i- 
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then Ay will scale as Aa^/'^; otherwise, Ay ^ Aa, and y{g,a) will be analytic on the 
G-curve. Taken with the condition (Jl^), the equation (]1^) has a unique solution: 



= ~ 0.0589, a, = 1 (19) 

V 288 

This will be seen to be the bicritical point anticipated above. At all the other points 
((7, ao((7)) on the G-curve (pTT]), y will be analytic, while F will be singular; we shall return 
to this point below. 

Now, the T-curve. Consider a point {g,a) at which alii{g,y) — y is analytic (the 
singular points have been treated in the preceding paragraph). By the implicit function 
theorem we know that y{g,a), as a root of the tree equation, will be analytic, unless 
the partial derivative of the equation with respect to y vanishes. The following equation, 
together with (P), defines the T-curve: 

aU2ig,yig,a))-l = (20) 

We solve the transcendental system (^ and (^) numerically. A convenient way to do so 
is to pick some value of g, solve Ili{g^ y)/Il2{g: y) — y = 0, and plug the value(s) of y into 
(^ or ( PPD to find a. In addition, (^) and ( PP] ) can be developed in powers of g near g = 0; 
there are two branches of the G-curve, a±{g). After a little algebra, we find 

63 

y±ig) = ±l-6gT—9^ + Oig^) 

1 3 9 (21) 
a,ig)^±---T-^9 + Oig^) 

Having canceled the first derivative of the tree equation to find the critical T-curve, we 
ought to search for a multicritical point on that curve that would be the result of canceling 
also the second derivative. However, we do not find this type of multicritical point realized 
in our system. 

This is the phase diagram of the TSMF model: 
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g 



As our system is symmetric under g — > — (7, we only show the right half. For complete- 
ness we show the a < quadrant, although only a > corresponds to real temperature. 
The bicritical point ((7^, a*) is marked. It is not only the solution of (|T^), but is also where 
the upper branch of the T-curve, aj^{g) meets the G-curve, aQ^g). For g > g*, there may 
well be solutions a+{g) above ao(fi'), but we cannot find them, as our generating function 
Il{g, A) breaks down in that regime: matrix models do not make sense above the geomet- 
ric transition. a_, on the other hand, continues indefinitely, approaching —1 as g ^ 00. 
So does the G-curve, although they never meet; in fact, their difference ao(fi') — a-{g) is 
asymptotically l/9g. 

Let us consider in more detail what happens to the solutions of the tree equation 
in the different regions of the phase diagram. First, fix some g > g*- For a < a-{g) or 
a > ao{g) the equation has no solutions: our solution for the free energy is not valid outside 
this region. For a_ < a < ao, y((7, a) is analytic. Near a_, the number of baby universes 
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on the tree diverges; y — y- ^ — a_, because we have canceled one derivative in y of 
the tree equation. Of the non-integral powers in F{g,a), the lowest, — a_, is canceled 
by the tree equation (^). Therefore F{g,a) ~ analytic + (a — a_)"^/^ + • ■ and a = \. 
The trees diverge as if the baby universes were not there, as | is the generic exponent 
for tree divergence, found, for instance, in the "Cayley" model (which just counts labeled 
trees with nothing at the nodes, in other words with Pn{g) = 1), or in one phase of the toy 
model of Ambj0rn, Durhuus, Frohlich and Orland (in which the role of a is played by 
the cosmological constant, and therefore 7str = \)- Near ao, the area of each baby universe 
diverges. Since (|T^) is not satisfied, y is analytic at However, since T{g, y^) = go, this 
is a critical point for F, which, according to (|TBp, scales as (ao — a)^/^ + analytic terms. 
Therefore 7str = the entire tree of baby universes diverges as if it were one, empty 
surface. 

Now fix some g < g*- For a > a_|_ or a < a_, the tree equation again has no solutions. 
At a± the tree diverges, and a = ^ as before. When a is between a_|_ and oq, there are 
two solutions to the tree equation; we must take the one which maximizes the free energy. 
By definition, when a = ao, one of these roots is yo, which satisfies r((7,|/o) = go and so 
would make the free energy diverge. However, this root is a minimum of the free energy, a 
"metastable" state not accessible to the system. Therefore, when g < g^, the area of the 
baby universes never diverges. (This is why this part of the G-curve is drawn as a broken 
line.) 

Finally, consider g = g*. As usual, a = | at a = a_. At a = ao = a-|-, the number 
of baby universes and their area diverge simultaneously, y is no longer analytic; as shown 
above. Ay scales as Aa^/"^. The lowest non-integral powers in the free energy, Aa^/"^ 
and Aa^/"^, are canceled by the critical conditions (p!6|) and (^Sj), respectively. Therefore 
F{g, a) ^ analytic + Aa^/"^, and 7str = 

Once again, we stress that there is no magnetization in the TSMF model. The T-curve 
is a tree-growing, not a spin-disordering, phase transition. The only exception is the a = 
line. Here, the a ^ and A ^ oo limits do not commute. If we take the a ^ limit 
first, the system will be magnetized for any value of A; taking the limits in reverse order, 
in other words putting a = in the TSMF model, yields no magnetization. 

One could repeat this calculation with a different matrix model potential U{(p) than 
(^. As long as the potential is still at the m = 2 point 0], the critical behavior of 

^ Strictly speaking it cannot be, since it is not even defined for a > ao; wliat is meant is that 
the series expansion of y for a < ao in powers of ao — a includes only non- negative integer powers. 
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the model should be qualitatively the same as that of the cubic model, with the same 
exponents, including 7str = ^ at the bicritical point. However, if the potential is tuned to 
one of Kazakov's m > 3 multicritical points, on the T-curve the generic a will still be |, 
but on the G-curve the generic 7str will now be — ^. At the bicritical point. Ay Aa™+i , 
and 7str = 

To conclude, we compare our results to related work. A number of years ago, several 
groups studied random surfaces embedded in continuous d-dimensional euclidean space 
(see and 0, and references therein). These models can be viewed as multiple gaussian 
systems coupled to a matrix model. It was argued that in the limit (i ^ oo, the surfaces 
are dominated by branched polymers, that is long thin tubes connected into trees. This 
situation can be viewed as a special case of the TSMF theory, its regime of small baby 
universes. This is not surprising, as the multiple gaussian models have only one coupling 
constant — the cosmological constant — which is coupled to the total size of the polymer, 
a role played by the matter temperature in the TSMF model. The reader will not be 
surprised to learn, therefore, that in certain regimes of the the multiple gaussian systems, 
the authors of [Qand 0find 7str = ^5 ^ result analogous, then, to our a = \. 

Recently, Durhuus and Ambj0rn obsereved that if a random surface model has 
7str = underlying, unbranched form, its "polymerized" version will have 7str = 

. This is in perfect agreement with our results, with m the multicritical index. This is 
because the trees of empty baby universes, which we found in |^ to dominate the A ^ oo 
limit, correspond precisely to the polymerized surfaces that Durhuus and Ambj0rn treat. 
There is another model, the (Tr(/)^)^ model introduced in ||^, which has phases with 
just the above exponents, and which Ambj0rn points out as an explicit realization of the 
construction in ||^. This is because the nonlocal term induces surface touching, giving 
surfaces that are trees of baby universes.^ The coefficient of the nonlocal term is therefore 
the analogue of the matter temperature in the TSMF model, and the phase diagrams of the 
two models are quite close, our bicritical point corresponding to Korchemsky's intermediate 
phase [ffO . 
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Appendix 

Here we give tiie expansion coefficients for tlie generating function Il{g, A), required 
in equation (jl^) and tlie subsequent anaiysis. From 0] we iearn tliat 

N a 3a^ + 6a + 2 1 , „ , 

U(g) = + - log 1 + 2a) 

3 (1 + (t) (1 + 2(7)2 2 ^ 

1/3 

T{g)= (3 (v/3(l- 34992(74)- 32 V)) 

Near tlie critical point qq = 1/ \l 108-\/3, ^[g) has the expansion (A^f = go — g) 
U{g) = TTo + 7riA(7 + n2Ag^ + n^/2Ag^/^ + 0{Ag^) 

I44.3I/4 

(3 + \/3)3 (2) 
_ 216(187+ 108^3) 

459 + 265v^ 
1152^2 • 3^/8(627 + 362^3) 



^3 1 , 



TTi = 



7r2 



^^^^ 5(362 + 209^3) 

We are interested in the expansion of Il{r{g, y)) about y = yoig); letting Ay = yo{g) — y 
and 'y = g/go^we get 

U{T{g, y)) = ^0 + ^lAy + ^aAy^ + ^^/^Ay^'"^ + ^(Ay^) 

TTo = TTo 

= %o7"^^^7ri (3) 

^2 = -^^7o7"'/'^i + 81(7^7-2/37r2 
7^5/2 = 243(7o'7"'/'7r5/2 

The rest of the terms in n((7, A) are of course analytic at y = yo{g), and their expansion is 
(j)Q + ^1 Ay + 4>2Ay'^ + (9(A?/^), where the coefficients are: 

-1/1 1 1 72/3\ 1 



^ 1__ 1 72/3\ 1 

2 V864^ ~ 14472/3 ^ 108 ~ 288 y ~ 3 



(70^ V8647 



V127 67V3 12 y 7I/3 

1 ^ 18^7o' 1 



272/3 ^2/3 2 
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The functions u{g), v{g), and w{g) can now be read off: 

u{g) = 0i(^) -7fi(^) 

v(g) = 2(M9) - Mg)) (5) 
5 

= -^^5/2(9) 
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